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^ ' Volume fluctuations are introduced in a statistical modelling of relativistic particle collisions. The 
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micro-canonical ensemble is used, and the volume fluctuations are assumed to have the specific scaling 
properties. This leads to the KNO scaling of the particle multiplicity distributions as measured in p+p 
interactions. A striking prediction of the model is a power law form of the single particle momentum 
spectrum at high momenta. Moreover, the mean multiplicity of heavy particles also decreases as a 



^ . function of the particle mass according to a power law. Finally, it is shown that the dependence 

of the momentum spectrum on the particle mass and momentum reduces to the dependence on the 

m 

^ ' particle energy. These results resemble the properties of particle production in collisions of high energy 

in 

C^^ ■ particles. 
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I. INTRODUCTION 

In collisions at relativistic energies many new particles are produced. Their number, masses 
and charges as well as their momenta vary from event to event. Most of the experimental 
results concern single particle production properties averaged over many interactions. It is well 
established that some of these properties, namely, mean particle multiplicities and transverse 
momentum spectra, follow simple rules of statistical mechanics. In proton-proton (p+p) colli- 
sions the single particle momentum distribution has an approximately Boltzmann form |lt] in 
the local rest frame of produced matter: 

where T, p and m are the temperature parameter, the particle momentum and its mass, re- 
spectively. At large momentum, ■p'^ ra, Eq. ([T]) gives: 

dN ( V\ / N 

Integration of ([T]) over momentum yields the mean particle multiplicity, (N), which is also 
governed by the Boltzmann factor for m» T: 

{N) r. ^mTfl- exp (- ^) . (3) 

The approximate validity of the exponential distributions dUlHl) is confirmed by numerous exper- 
imental results on bulk particle production in high energy collisions. The agreement is limited 
to the low transverse momentum [px < 2 GeV) and the low mass {m < 2 GeV) domains. 
However, the temperature parameter T extracted from the data on p+p interactions is in the 
range 160-190 MeV [2]. Thus, almost all particles are produced at low pt and with low masses. 

Along with evident successes there are obvious problems of the statistical approach. The 



probability 



Refs. 
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P{N) to create N particles in p+p collisions obeys the KNO scaling [2,| (see also 



]), namely: 

P{N) = {N) -1 i;{z) , (4) 

where (A^) is the mean multiplicity and the KNO scaling function tp{z) only depends on z = 
N/{N). The mean multiplicity increases with increasing collision energy, whereas the KNO 
scaling function remains unchanged. The latter implies that the scaled variance uj of the 
multiplicity distribution P{N) grows linearly with the mean multiplicity: 

^ = ^ \^^^ ^ oc (AT) . (5) 
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A qualitatively different behavior is predicted within the existing statistical models [Tl, [sl, Isj, 
namely the scaled variance is expected to be independent of the mean multiplicity: 

u) ~ const ^ 1 . (6) 

This contradiction between the data and the statistical models constitutes the first problem 
which will be considered in this paper. 

The second and the third problems which will be addressed here concern particle production 
at high (transverse) momenta and with high masses, respectivel y. I n these regions the single 



particle energy distribution seems to obey a power law behavior 10|] : 



At p ^ m, Eq. (JTj) gives: 



'" P-"' . (8) 



with Kp = K. Integration of ([7]) over particle momentum yields the mean multiplicity which 
follows a power law dependence on the particle mass: 

(A^) ~ m"^™ , (9) 

with Km = Kp — 3. The Eqs. ([8]) and ([9]) approximately describe the data on spectra of light 
particles at large (p > 3 GeV) (transverse) momenta and on the mean multiplicity of heavy 
(m > 3 GeV) particles, respectively. The power law parameters fitted to the data are Kp = 8 
and Km ~ 5 [10]. One observes a growing disagreement between the exponential behavior ([2]), 
([3]) and power law dependence (E]), (E]) with increasing (transverse) momentum and mass. At 
p = 10 GeV or m = 10 GeV the statistical models underestimate the data by more than 10 
orders of magnitude. 

In the present paper we make an attempt to extend the statistical model to the hard domain 
of high transverse momenta and /or high hadron masses (hard domain). The proposal is inspired 



by statistical type regularities 



lOJ l in the high transverse mass region, as well as by the recent 



work on the statistical ensembles with fluctuating extensive quantities HI]. We postulate that 
the volume of the system created in p+p collision changes from event to event^. The main 
assumptions of the proposed approach are the following. 



^ The volume fluctuations in hadron statistical physics were first introduced in the framework of the isobaric 
ensemble in Refs. [12|, |l3|. 



1. Each final state created in p+p interactions is identified with a micro-state of a micro- 
canonical ensemble (MCE) defined by the volume V, energy E, and conserved charges 
of the system. By definition of the MCE, all its micro-states appear with the same 
probability. 

2. The volume of micro-canonical ensembles represented in p+p interactions fluctuates from 
collision to collision. The volume probability density function, Pa{V), is given by the 
scaling function, PaiV) = V (paiV/V), where V is the scaling parameter. 

The model based on these assumptions will be referred as the Micro-Canonical Ensemble with 
scaling Volume Fluctuations, the MCE/sVF. Calculations presented in this work aim to 
illustrate the main idea and therefore they are performed within the simplest possible model 
which preserves the essential features: an ideal gas of massless particles and a mixture of 
massless neutral with heavy neutral particles. The Boltzmann approximation has been used. 

The paper is organized as follows. The description of the relativistic gas of particles within 
the MCE is presented in Section II. The scaling volume fluctuations within the MCE framework 
are introduced in Section III. Then the basic equations of the MCE/sVF are given. Next, the 
results on the multiplicity distribution, the single particle momentum spectrum and the mean 
multiplicity of heavy particles are obtained. The paper ends with the summary and closing 
remarks presented in Section IV. 

II. MICRO-CANONICAL ENSEMBLE 

The MCE partition function for the system with N Boltzmann massless neutral particles 
reads U,\l^: 

where E and V are the system energy and volume, respectively, g is the degeneracy factor, 
and A = gVE^/ir'^. The MCE partition function fllUI) includes exact energy conservation, but 
neglects the momentum conservation. The MCE multiplicity distribution is given by: 

P (N-EV) = ^"^^'^^ (11) 



W{E, V) is the total MCE partition function 

W{E,V) ^ X;^^(^'^) = ^ °^^('3'3'2'^j ' (12) 



N=l 



where o-^s is the generahzed hyper-geometric function (see Appendix). For A ^ 1 the mean 
multiphcity equals to [8|]: 



(N). 



Y,N Pmce{N-E,V) = {A/27) 



1/4 



(13) 



N=l 



where Pmce{N] E, V) was approximated by the normal distribution [8]: 



P^ce{N;E,V) ^ (27r^„,e-(iV) 



mce ) exp 



( N - {N)mce y 



2uJr 



{N)r 



(14) 



with Umce = {{N^)mce — {N)mce)/{^)mce = 1/4. Note that in the grand canonical ensemble 
(GCE) the multiplicity distribution is equal to the Poisson one: 



Pg,e{N;T,V) 



N 



N\ 



exp(-A^) 



(15) 



where A^ is the mean multiplicity in the GCE. The distribution 01 5p approaches the Gaussian 
for large A^: 



PgUN;T,V) ^ {2'Kujg,,-NY^" exp 



[N - N) 

2 id gee ■ N 



2 1 



(16) 



with Ugce = {N^-N )/N = 1. 

The numerical calculations presented in this paper will be performed for g = 1 and the 
energy density which corresponds to the temperature parameter T = 160 MeV. The latter 
relates the values of E and V via equation: 
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E 



V T" . 



TC^ 



(17) 



The mean multiplicity {N)mce in the MCE ( [T3|) is then approximately equal to the GCE value: 

N = ^V T"^ . (18) 

The approximation {N)mce = N is valid for A^ ^ 1 and reflects the thermodynamic equivalence 
of the MCE and the GCE. The scaled variance of the MCE distribution is Umce = 1/4 8J], and is 
approximately independent of {N)mce already for {N)mce > 5. Thus, despite of thermodynamic 
equivalence of the MCE and GCE the value oiumce is four times smaller than the scaled variance 
of the GCE (Poisson) distribution, Ugce = 1- 

The single particle momentum spectrum in the GCE reads: 



Fgceip) = 



1 dN 



V 



N p^dp 27r2 N 



exp 



p_ 
T 



2T3 



exp 



p_ 
T 



(19) 
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FIG. 1: (Color online) (a): The multiplicity distribution of massless neutral particles in the MCE 
(jlip . dashed line, and the GCE (jlSp . dotted line, (b): The momentum spectrum of massless neutral 
particles calculated within the MCE (|2Up . dashed line, and the GCE (J19p . dotted line. The system 
energy is ^ = 3iVT = 24 GeV for both plots. 



whereas the corresponding spectrum in the MCE is given by: 



" race \P) 



dN 



V 



{N)mce P^dp 27r2(iV) 

1 

2^3 



E 

N=2 



Wn-i{E-p,V) 



V 



W{E,V) 



2tt^N), 



{N)r 



1 ^ A^(3iV-l) 

N=2 



(3iV-4)! 



1 - 



p_ 

E 



3N-4 



Pmce{N;E,V) 



f{p;E,V) 

(20) 



where f{p;E,V) is the MCE analogue of the Boltzmann factor, exp{—p/T). From Eq. ( 120|) 
follows f{p = 0;E,V) = l. Both (IT9l) and (!20l) are normalized such that J^ p^dp Fgce{p) = 1 
and /q p'^dp Fmcdp) = 1- 

Figure [T^ shows a comparison of the MCE and GCE results for the multiplicity distribution 
and Fig. [T)d shows the momentum spectrum for A^ = 50. The MCE spectrum is close to 
the Boltzmann distribution ( IT9|) at low momenta. This can be shown analytically using the 
asymptotic form of the generalized hyper-geometric function (see Eq. (1A4p in Appendix) at 
E ^ oo and p/E < 1: 



f{p;E,v) 



W{E-p,V) _ (E-p) 

W{E,V) 



2oi^3(; ^A.2 



E^ 



0F3 ( 



.459. V{E-pf 
' 3' 3' ' 277r2 



.459. VE3 \ 
' 3' 3' ' 27iT^) 



exp 



p_ 
T 



(21) 



The MCE spectrum decreases faster than the GCE one at high momenta. Close to the threshold 
momentum, p = E, where the MCE spectrum goes to zero, large deviations from (TT^ are 



observed. In order to demonstrate this the MCE and GCE momentum spectra are shown in 
Fig. [lb over 90 orders of magnitude. 

III. MCE WITH SCALING VOLUME FLUCTUATIONS 

Let us consider a set of micro-canonical ensembles with the same energy E but different 
volumes V. The probability density which describes the volume fluctuations is denoted by 
Paiy). The distribution of any quantity X can then be calculated as: 



/•oo 
P^{X-E) = / dV P^{V)Pmce{X-E,V) , 

Jo 



(22) 



where PmceiX;E,V) is the distribution of the quantity X in the MCE with fixed E and V. 
Further more it is assumed that the distribution Pa{V) has the scaling form: 

Pa{V) = V~' MV/V) , (23) 

with the volume parameter V being the scale parameter of Po,{V). From Eq. fpHl) follows that 
the scale parameter is proportional to the average volume and the proportionality factor follows 
from the normalization conditions of the scaling function (see below). Equations ( 122|) and (123|) 
are the basic equations of the MCE/sVF. 

Let us introduce an auxiliary variable y defined as y = (V/V) . Using Eqs. (fT3|l and 
flTTl \TE\i the mean multiplicity of massless particles in the MCE (1131) can be written as: 

{N)mce = Ny, (24) 

where 

The average temperature T is related (see Eq. (ITTl) ) to the energy and volume parameter V as: 

E = ^VT\ (26) 

Within this paper T = 160 MeV is assumed. The change of the volume aX E = const leads to 
the corresponding change of the system temperature: 

T = T-(y/Vf'^ = T/y. (27) 

The temperature fluctuations resulting from the volume fluctuations will be discussed below. 



Using the variable y, the volume integral in Eq. fl22l) can be conveniently rewritten as: 

P^{X-E) = {yy^ / dV<P^ {V/V) Prr,UX;E,V) = / dy My)F^UX;E,yW) , (28) 
Jo Jo 

where ipa{y) = '^y^4'a{y^)- Choosing ipct{y) = S{y — 1) one recovers the MCE with V = V and 

\-' ' /mce -' ' • 

The scahng function ipaiv) will be required to satisfy two normalization conditions: 



dy i/Ja{y) = 1, / dyyipaiy) = 1. (29) 

Jo 

The first condition guarantees the proper normalization of the volume probability density func- 
tion, Jg dVPa{V) = 1. The second condition is selected in order to keep the mean multiplicity 
in the MCE/sVF equal to the MCE mean multiplicity. 

A. KNO-scaling 

The multiplicity distribution in the MCE/sVF is: 

/•OO /"OO 

PaiN;N) = / dV PaiV) P^UN;E,V) = / dy Mv) Pn^ceiN; E,y^V) . (30) 

Jo Jo 

At iV :^ 1, the particle number distribution flHUl) can be approximated as: 

Pa{N-N) = {N)~'Mz) , (31) 

where z = N/{N)a , and the mean multiplicity {N)^ is given by: 

(iV), = J2 N Pa{N;N) = N / dyytPM = N ■ (32) 

The approximate equality of {N)a and A^ is satisfied for A^ ^ 1 due to the second normalization 
condition ( l29l) . The KNO scaling of the multiplicity distribution Pa{N;N) follows from the 
assumption of the scaling of the volume fluctuations fl2^ . 

For convenience, a simple analytical form of the scaling function, ipa will be used: 

My) = ^y'~' ^M-ky), (33) 



where T{k) is the Euler gamma function. It was found 13| that the function flHHj) with k = 
4 approximately describes the experimental data on KNO scaling in p+p interactions. All 
numerical calculations obtained within the MCE/sVF and presented in this paper will be done 



using the function ( l33l) with k = 4. Note, that the function (1331) satisfles both normalization 

8 



conditions (!29l) for any A; > 0. In order to check the sensitivity of the results derived within 



the MCE/sVF to the shape of the KNO scahng function, the calculations presented 
repeated using the scaling function resulting from the fit to the experimental data 

i'aiy) = ay" exp(-6y2) , 



Delow were 

4|: 



(34) 



where the values of the parameters are the following: a = 1.19, 6 = 0.62 and c = 0.66. 
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FIG. 2: (Color online) (a): The scaling functions determining the volume fluctuations discussed in this 
paper. Solid line indicates the function (I33p with k = 4 for which all presented results are calculated, 
whereas the dashed line shows the function (I34p used for the sensitivity tests, (b): The multiplicity 
distributions calculated within the MCE/sVF for three values of mean multiplicity. The calculations 
are performed for the function (jSSp with k = 4. The distributions are presented in the KNO variables, 
{N) ■ PiN) and z = N/{N). 

In Fig. [2^ both KNO scaling functions (!33|) and (!M|) are plotted for a comparison. The qual- 
itative behavior is similar, however there are substantial quantitative differences. Nevertheless, 
the main results presented below are the same for both scaling functions. 

The multiplicity distributions calculated within the MCE/sVF for the three mean multiplic- 
ities, A^ = 5, 10 and 50, are shown in Fig. [2b in the KNO variables. The scaling function ( l33l) 
is also shown for a comparison. These numerical results indicate that even for relatively low 
mean multiplicities the KNO scaling is approximately obeyed by the MCE/sVF distributions. 

The KNO scaling of the multiplicity distribution implies that the scaled variance of the 



distribution increases in proportion to the mean multiplicity. For A^ ^ 1 one gets: 



Un 



(N). 



K {N)^ , 



(35) 



where k = const > 0. Equation (1351) results from an approximation of {N'^)a for A^ ^ 1, 
namely: 

CO «QQ 

{N')^ = J2N'P^iN;N) = n' / dy y' Mv) = (l + ^^) {N)l . (36) 

The positive value of k in Eq. ( l36l) follows from the normalization conditions (l29l) : 

1 /»oo 

dyy^My)-'^= / dyiy-lfMy) > 0, (37) 

For the function ^ipaiy) defined by Eq. ( l33l) one finds, k = k~^. 
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FIG. 3: (Color online) A comparison of the multiplicity distributions of massless neutral particles 
calculated with the MCE/sVF (solid line) and the MCE (dashed line). The system energy is £' = 
3lVr = 24 GeV. 

In Fig. [3] the multiplicity distributions obtained within the MCE and the MCE/sVF for 
A^ = 50 are compared. The scaled variance of the MCE/sVF distribution for A^ = 50 is about 
12.5, whereas the scaled variance of the MCE distribution is 1/4. This large difference in the 
width of the MCE/sVF and the MCE distributions is clearly seen in the figure. 

The volume averaging procedure results in the mean value given by: 



{V)a 



dVVP^iV) = V / dyy^My) = aV . 
Jo 



(3^ 
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For ipa (1331) one gets: 

{k + l){k + 2){k + 3) 



(39) 



fc3 

which gives a = 3.28 for k = 4. Thus, in general, the mean volume in the MCE/sVF is not 
equal to the corresponding MCE volume. 

B. Power Law in Momentum Spectrum 

The volume fluctuations in the MCE/sVF significantly increase the width of the multiplicity 
distribution. They are also expected to modify the single particle momentum spectrum. This 
is because for a fixed system energy, the volume of the system determines the energy density, 
and consequently, the effective temperature of particles. 

The single particle momentum spectrum within the MCE/sVF can be directly calculated 

from Eq. fl22l) and it reads: 

1 / dN \ 1 r°° V 

Faip) = J^i^ = 7^ / dVP^{V) — f{p;E,V) 

{N}a \p^dp/a {N}a Jq 2tX^ 

= -A -^y^ -^7^^ f^ (1 - ^) PJN-N). (40) 

{N)^ 2^3^^ (3A^-4)! \ e) "V , ; v ; 



The formal structure of the expression ( 1401) is similar to the structure of the corresponding 
expression derived within the MCE (1201) . The only, but the crucial, difference is that the narrow 
MCE multiplicity distribution used for averaging the particle spectrum in Eq. (1201) is replaced 
by the broad MCE/sVF multiplicity distribution in Eq. (HUj) . The spectrum Fa{p) fulfills the 
normalization condition, J^ p'^dp Fa{p) = 1. From Eq. (HUj) one finds Fa{p = 0) = a-Fg^eip = 0). 
Equation ([39]) gives a = 3.28 for the tpa function (ESD with k = 4. Thus, in the MCE/sVF 
there is an enhancement of the momentum spectrum at p — ^ compared to the GCE and MCE 
results. 

The single particle momentum spectrum calculated with the MCE/sVF for the volume 
scaling function flHHl) is shown in Fig. H^. A striking new feature of this spectrum is the 
presence of a long power law tail. In the momentum range from several GeV to about 20 GeV 
the spectrum can be approximated by: 

Fa{p) = Cpp-""" , (41) 

with Cp and Kp being the normalization and power parameters, respectively. For momenta 
smaller than 3 GeV the spectrum starts to deviate significantly from the power law parametriza- 
tion and its local inverse slope parameter is close to the temperature of the corresponding GCE, 
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FIG. 4: (Color online) (a): The momentum spectrum of massless neutral particles calculated within 
the MCE/sVF (gO]), soHd Hne, and the GCE ([19]), dotted hne. The approximation (gSD of the 
MCE/sVF spectrum is shown by the dashed-dotted line, (b): The same as in the left panel but 
only for the low momentum region. The approximation (I42p is not shown as it overlaps with the 
MCE/sVF line. The system energy is E = 3NT = 24 GeV for both plots. 

T = 160 MeV. This is shown in Fig. HJd, where the momentum spectrum for p < 3 GeV is pre- 
sented. A rapid decrease of the spectrum starts at p > 20 GeV, when the threshold value 
p = 24 GeV is approached. The above features of the MCE/sVF momentum spectrum resem- 
ble features of the transverse momentum spectrum of hadrons produced in high energy p+p 
interactions [Iq . 

The power law dependence fHT]) of the momentum spectrum at high momenta can be derived 
analytically, namely: 

V 



Fa{p) 



Fmce{p)i'a{y)dy 



2'K^N 



*'r(* + 4)_.«|^^_^^_,_, ^ 



/ dyipa{y)y^ exp (- ^?/ 
11.27 GeV^ (p + 4T)-^ 



(42) 



2r(fc) 

where T = 160 MeV and k = 4 is set in the last expression. Note that according to Eq. fl7r|) 
the temperature, T = T/y, fluctuates^ in the MCE/sVF due to the volume fluctuations. The 
volume fluctuations are responsible for the appearance of the power law tail. The temperature 



^ In Refs. Il6| it was shown for the first time that the properly chosen temperature fluctuations may lead to 
the Tsallis distribution [ITil with power law spectrum at high (transverse) momenta, see also Ref. il^] . 
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probability distribution can be easily derived: 



P(r) = |;^WT/T) = ^^ 



k+1 



exp 



4 

T 



The result is shown in Fig. [5]d. 
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FIG. 5: (Color online) (a): The momentum spectrum of massless neutral particles calculated within 
the MCE/sVF ()40p for three values of the mean multiplicity equal to A^ = 25, dashed-dotted line, 
A^ = 50, dashed line, and N = 100, solid line. The system energy is E = 3NT. (b): The temperature 
probability distribution (03]) with k = 4 in the MCE/sVF. 



The analytical approximation in Eq. fH2l) is valid for p <^ E. It describes the particle 
spectra at p <^ E including the low momentum region (see Fig. HI (a)). For p ^ T one gets 
the power law parameters ( 14T1) : Cp = 11.27 GeV^^ and Kp = 8. The power law dependence of 
the particle energy spectrum appears when there is a distinct region of the particle energies in 
which T <^p <^ E. 

Equation fH21) leads to an exponential spectrum for p <^T : 



2T 



To 



(44) 



with the inverse slope parameter Tq = T/2 and the normalization constant a = 3.28. For 
arbitrary k > 0, the parameter a is given by Eq. (l39ll . and Tq = Tk/{k + 4). Within the 
used approximation the power parameter, Kp, depends only on the form of the volume scaling 
function, whereas the normalization constant Cp in Eq. (HTj) is sensitive, in addition, to the 
temperature T. 
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Note that the scahng volume fluctuations at a flxed system energy lead alone to the KNO 
scaling and the power law behavior of the spectra. Systems with small volume, V -^ V, 
have predominately a small particle multiplicity, N <^ N. Due to the energy conservation 
these systems have large temperatures, T ^ T, which results in the power law behavior of 
the momentum spectra and the mean heavy particle multiplicity (see the next Section). The 
energy fluctuations at flxed volume do not lead to the power law dependence. This is because, 
systems with a small energy, E ^ E, have predominately a small particle multiplicity, N <^ N, 
and consequently the temperature is approximately independent of the system energy, volume 
or multiplicity. Thus, no power law behavior appears in this case. 

On the other hand, the energy fluctuations can be added into the present consideration. 
They do not change the principal behavior of particle spectra. In order to illustrate the in- 
dependence of the power law tail of the system energy the MCE/sVF spectra calculated for 
three different values of E and compared in Fig. [5ti. In the overlapping power law region (left 
to the cuts from the energy conservation) the spectra are similar, i.e. they do not depend 
on the average multiplicity. This implies that averaging over collisions with different energies 
deposited for particle production would not change the particle spectra for momenta below 
the region affected by the steep threshold decrease. In particular, it follows from Fig. [5H that 
the energy fluctuations in the region ii^ = 10 ^ 50 GeV do not change the momentum spectra 
at p < 10 GeV. The same arguments are applied to the mean multiplicity of heavy particles 
discussed in the next Section and shown in Fig. [6t) 

Finally, it is important to note, that the momentum spectrum for particles with non-zero 
mass, m, obeys a power law dependence in the particle energy. This is obvious from Eq. fH2|) . 
because the particle momentum, p, in the exponent can be replaced by the particle energy. 



\/p'^ + m?, resulting in: 

^ r^ / . ^ 4 I v^ 



-^"(^) - ':r^Tf / dyi)a{y)y exp 

zvr^A' Jo 



m? 



T 



= ^ p.^x T {^p^ + m^ + Tk) \ (45) 

Consequently, the dependence of the momentum spectrum on the particle momentum and 
mass reduces to the dependence on the particle energy. This energy scaling predicted by the 
MCE/sVF seems to coincide with the scaling in the transverse mass of the quarkonia spectra 



measured in p+p interactions [10 1. 
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C. Power Law in Mass Dependence 

The volume fluctuations together with the specific scaling function fl33l) lead to the power 
law dependence of the momentum spectrum at high momenta. It may be expected that these 
volume fluctuations will also cause a similar dependence of mean multiplicity of heavy particles 
on the particle mass. This subject is discussed below in the approximation of heavy and thus 
non-relativistic particles. 

The MCE partition function for n heavy particles is given by |l4l . [l9| : 

1 / GV \" m^ ,^ x32i_i , , 



where m is the mass of the heavy particle and G is its degeneracy factor. Using Eqs. ( TTOl) and 
P6l) one can calculate the MCE partition function for A^ massless and n heavy particles. It 
reads pj,|l9|: 

ZN,niE, V) = dE, dE2 W^iE,, V) nrXE2, V)5[E - E, - E2] 

Jo Jo 

1 1 f qV\ ( GV \" m^ ,^ xqAT-uS" 1 / _N 

In what follows the values of the degeneracy factors are assumed to he g = G = 1. Note, 
that due to the exact energy conservation, the MCE partition functions, Wn{E, V) flTUj) and 
^n{E, V) ( l46l) . are defined for a particle number larger than zero, A^ > 1 and n > 1. However, 
the MCE partition function, Z^^niE, V) f HTj) . requires only N + n > 1, so that either A^ or 
n can be equal to zero. The mean heavy particle multiplicity in the system with two types of 
particles, massless and heavy, is: 



^ 00 



ZiE, V) 

^ ' ' N,n=0 

where Zn^E. V) is given by Eq. gTl), Z{E, V) = E",n Zn^E, V), and Zo,o(^, V) = 0. The 
energy conservation implies that there is a maximum number of heavy particles, rimax = [^] • 
Figure EK shows the dependence of the mean multiplicity of heavy particles on mass m calculated 
within the MCE fHHj) . The plot starts at m = 1 GeV. This is because at lower masses, 
m = T = 160 MeV, and thus the non-relativistic approximation used for heavy particles is not 
justified. Furthermore, if m < 1 GeV, a contribution of heavy particles to the system energy 



density can not be neglected, and this would change the temperature parameter T 19|. An 
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exponential decrease with the inverse slope parameter of about 160 MeV observed for masses 
of several GeV, is followed by a steepening decrease towards the threshold. 

At low masses, m <^ E, the MCE mean multiplicity is close to the corresponding multiplicity 
calculated within the GCE, see Fig. (a): 



n 



\/ (mT/27r)3/2 ^^^ f ^^ _ ^^g^ 



This can be shown analytically as follows. For the considered values of the total energy E 
and volume V the average multiplicity {n)mce is essentially smaller than one. Thus, the main 
contribution to the numerator of Eq. (HHl) comes from ZN^n=i, while the main contribution to 
the denominator comes from ZN^n=o, thus: 

^ , , T? f 1 5 7 V(E-m)^\ 

,. ^ En=o Zn,i _ V2my^VE^r^ °^H' ^' ^^ ^^ ^^^j .,n^ 

\ri'/mce - ^oo ^ - p2 p / 4 5 o V E'^\ ' ^"^^^ 



Using Eq. flA4p from Appendix one can show that in the region T <C m <C -E it follows: 

{n)mce = n. (51) 

Note that Eq. fl3Ul) is valid not only for {n)mce ^ 1, but for all values of {n)mce including also 
{n)mce ^1- In the limit VE^ -^ oo the distribution of heavy particles in the MCE has the 
Poisson form Pmceiri) = exp(-rr) n"/(n!), for nm -C E. Thus, {n)mce = -Pmce(l)/-P»nce(0) = n 
at all n <^ E/m. 

The mean multiplicity of heavy particles within the MCE/sVF is given by: 



{n)a = / dy ^paiy) {n)mce , (52) 

and its dependence on the mass m is shown in Fig. [6] (b). This dependence resembles the 
behavior of the MCE/sVF momentum spectrum for massless particles shown in Fig. |H 
The mean multiplicity (n)^ as a function of m obeys the power law behavior: 

(n)^ = N Cn m-^™ , (53) 

in the large central interval of heavy particle masses. At high m, the (n)„ is larger than 
corresponding MCE and GCE values by many orders of magnitude. This shows non-equivalence 
of the MCE/sVF and the GCE or the MCE for mean multiplicity of heavy particles even in the 
thermodynamic limit. Another example of such non-equivalence was already obtained for the 
mean volume, see Eq. ( l38l) . For the masses close to the threshold mass, -m = 24 GeV, the mean 
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FIG. 6: (Color online) (a): The mean multiplicity of neutral heavy particles, (n), as a function of 
heavy particle mass, m. Note, that the horizontal scale starts at m = 1 GeV. The mass dependence 
obtained with the MCE (|15|) is indicated by the dashed line, whereas the GCE dependence P^ is 
shown by the dotted line. The system energy is E = 3NT = 24 GeV. (b): The mass dependence (I52p 
within the MCE/sVF is shown by the solid line, and the dependence for the GCE (I49p is indicated 
by the dotted line. The approximation (j54p of the MCE/sVF results is shown by the dashed-dotted 
line. The system energy is -E = 3A^ T = 24 GeV. 

multiplicity steeply decreases to zero. The power law approximation of the mass dependence 
can be derived analytically: 



{n)a 



dyipaiy) {n)r 



^X'vfdyMy) y'" 



N 



v^ k^ T{k + 5/2) T' 



;A;-fl 



m 



3/2 



A^ ■ 0.81 • GeV^ 



exp 

^3/2 



m y 



(54) 



2V2r{k) {m + Tkf+^i^ " """ "" {m + Tkyy^ 

In the last expression T = 160 MeV and k = 4 are set. The approximation in Eq. (l54|) is 
valid for the particle masses much smaller than the system energy, m <^ E. The analytical 
expression ( IMj) describes the numerical results on particle multiplicity at 1 GeV< m <^ E. This 
is illustrated in Fig. [Ub, where both are compared. If m ^ T, the power law dependence fl3^ 
of the mean multiplicity on m appears with Cm — 0.81 GeV^ and K^ = 5. This approximation 
is valid when there is a distinct region of m in which T <^m <^ E. 

The relation Km = Kp — 3, suggested by the experimental data and derived using a simple 
statistical considerations [lO], is valid within the MCE/sVF. Note, that also the values of the 
power law parameters, Kp = 8 and Km = 5, are close to the corresponding parameters ex- 
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tracted from the experimental data [10]. Within the MCE/sVF these parameters are, however, 
dependent on the choice of the volume scaling function. For the function ( IMj) one gets Kp = 5.7 
and K^ = 2.5. Thus, a quantitative comparison with the data requires further studies. 

IV. SUMMARY AND CLOSING REMARKS 

In this paper the statistical approach to particle production in high energy collisions is 
extended. A model called the Micro-Canonical Ensemble with scaling Volume Fluctuations, 
or the MCE/sVF, is proposed which incorporates volume fluctuations in the micro-canonical 
ensemble. A class of scaling volume fluctuations is considered. Analytical and numerical 
calculations are performed within the simplest formulation of the model which preserves its 
qualitative features. 

First, it is shown that the model leads to KNO scaling of the multiplicity distribution. The 
volume scaling function is selected in order to reproduce the KNO scaling function measured 
in p+p interactions. 

Second, the single particle momentum spectrum is calculated for massless particles. The 
MCE/sVF spectrum has features which resemble properties of the transverse momentum spec- 
trum of hadrons produced in very high energy collisions. In particular, a long power law tail 
appears as a result of the volume fluctuations. 

Third, the dependence of the mean multiplicity of heavy particles on the particle mass is 
obtained. The mass dependence approximately obeys a power law decrease, similar to what 
is observed experimentally. The power parameters for the momentum spectrum and the mass 
dependence of mean multiplicity are related as Km — Kp — 3, again, similar to experimental 
data. 

Fourth, it is also shown that for particles with arbitrary mass the dependence of the momen- 
tum spectrum on the particle mass and momentum reduces to the dependence on the particle 
energy. 

A quantitative comparison with the experimental data requires a significant additional work, 
in particular the introduction of the proper degrees of freedom and all related conservation laws. 

For more than 30 years the production of hadrons at high transverse momenta and/or with 
high masses, which follows the power law dependence, is discussed within dynamical QCD based 
models of parton scattering with a large momentum transfer. In particular, the perturbative 
QCD calculations quantitatively describe many experimental results. Thus, in view of the 
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evident success of the QCD, the apphcabihty range of the MCE/sVF and a relation of the 
model to the QCD based approaches should be critically studied. 

Clearly, further tests of the MCE/sVF are needed. Other qualitative features of p+p inter- 
actions should be discussed. The most important is the jet-like correlations of high transverse 
momentum particles. These correlations are usually described within the perturbative QCD 
for hard parton scattering accompanied with phenomenological parton-hadron fragmentation 
function. It is not yet clear whether the jet structure of particle production at high transverse 
momenta can be reproduced within the MCE/sVF. Clearly this would require an introduction 
of the momentum conservation. The multiplicity distributions in p+p collisions obey some 



additional regularities, e.g., the Golokhvastov scaling of the semi-inclusive spectra 20|. They 
should be also studied within the MCE/sVF. 

An effort to extend the model to nucleus-nucleus collisions is needed. First predictions can, 
however, be obtained from Eqs. (I42l) and (15^ . Namely, the particle momentum distribution at 
high (transverse) momenta and the mean multiplicity of heavy particles are expected to behave 
as: 

"^^ - ' C, (N) p-^^ , (55) 



p^dp 

(n) ^ Cm {N) m-^- , (56) 

where parameters Kp and K^ in Eqs. ( I55l 1561) depend only on the form of the volume scaling 
fluctuations, and the parameters Cp and Cm depend, in addition, on the mean temperature. 
Under the assumption that the volume scaling function and mean temperature are energy and 
system size independent, the MCE/sVF predicts the proportionality of dN/{p^dp) at large p 
and of (n) at large m to the average total multiplicity, (A^). Data on system size dependence of 
the mean multiplicity of J /if) mesons as well as the particle yield at high transverse momenta 



seem to confirm this prediction 2l|, |22[ |. 



Acknowledgment s 

We would like to thank F. Becattini, A. S. Golokhvastov, W. Greiner, M. Hauer, I. N. Mishus- 
tin, St. Mrowczynski and P. Seyboth for numerous discussions. This work was in part supported 
by the Program of Fundamental Researches of the Department of Physics and Astronomy of 
NAS, Ukraine and the Institute VI - 146 of Helmholtz Gemeinschaft, Germany. V.V.B. would 
like also to thank for the support of The International Association for the Promotion of Coop- 



19 



eration with Scientists from the New Independent states of the Former Soviet Union (INTAS), 
Ref. Nr. 06-1000014-6454. 

APPENDIX A 

The generahzed hyper-geometric function, also known as the Barnes extended hyper geo- 
metric functions defined by the following series [23;j ] : 

pFg{ai, 02, ... , ap] 61, 62, • • • , bq] z) = y , (Al) 

^ [Ol)k[02)k--- [Oq)k K\ 

where (a)^ is the Pochhammer symbol: 

(a)fe = ^^^^^ = a(a + l)...(a + fc-l) . (A2) 

r(a) 

The Euler gamma function T{x) has a simple form for integer k and halfinteger k + 1/2 argu- 
ments: 

r(fc) = (A:-l)!, r^fc + l) = ^'^'•••^'j^^'^^ v^- (A3) 

The asyrnptotic behavior at z ^ 00 of the hyper-geometric function pFq for p=0 and q=3 is 



given by |23|: 
.F3(;6:,i.6,;--) ^ ^^^^IMIMexp[4^..K|-'.— -0(1 + o(i=)). (A4) 
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